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^ ■ Abstract 

We employ the conformal bootstrap to re-examine the problem of finding 
the critical behavior of four-Fermion theory at its strong coupling fixed point. 
5_j ■ Existence of a solution of the bootstrap equations indicates self-consistency of 

the assumption that, in space-time dimensions less than four, the renormal- 
ization group flow of the coupling constant of a four-Fermion interaction has 
a nontrivial fixed point which is generally out of the perturbative regime. We 
exploit the hypothesis of conformal invariance at this fixed point to reduce 
the set of the Schwinger-Dyson bootstrap equations for four-Fermion theory to 
three equations which determine the scale dimension of the Fermion field ip, the 
scale dimension of the composite field ipip and the critical value of the Yukawa 
coupling constant. We solve the equations assuming this critical value to be 
small. We show that this solution recovers the fixed point for the four-fermion 
interaction with iV-component fermions in the limit of large N at (Euclidean) 
dimensions d between two and four. We perform a detailed analysis of the 
1/iV-expansion in d = 3 and demonstrate full agreement with the conformal 
bootstrap. We argue that this is a useful starting point for more sophisticated 
computations of the critical indices. 
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1 Introduction 



Conformal symmetry governs the critical behavior of a quantum field theory near the 
fixed points of its renormalization group flow ]I], ^ |3|. It is most powerful in two 
spacetime dimensions where the requirement of conformal invariance of correlation 
functions leads to a solution of a quantum field theory in terms of a few parameters ||. 
In this Paper we shall show that conformal symmetry is also useful in higher dimen- 
sions. In any spacetime dimensions, it determines the form of the two and three-point 
Green functions up to a few constants, the scaling dimension of the field operators 
and the value of the coupling constant at the fixed point. When this information 
is combined with Schwinger-Dyson equations, it can provide a powerful tool for the 
analysis of field theories. We shall apply the conformal bootstrap technique to four- 
Fermion interactions in dimensions 2 < d < 4. An application of this technique to 
A0 4 theory || and some preliminary work for spinors with four-Fermion interactions 
H between 2 and 4 dimensions have been given previously by Yu. Makeenko. 

Fermionic field theories with four-fermion interactions have a long history. A 
classic example is the Nambu-Jona-Lasinio model which was proposed some thirty 
years ago by Nambu and Jona-Lasinio |7| (see also ||) as a quantum field theory which 
exhibits dynamical chiral symmetry breaking for sufficiently strong attractive fermion- 
anti-fermion interactions in four (space-time) dimensions. The study of related models 
has recently seen a revival with speculations that the symmetry breaking mechanisms 
which it describes could be used in the standard model to describe the Higgs particle 
ClS db top quark condensate || |i~0f . 

In the conventional perturbative expansion in powers of the coupling constant, 
four-fermion interactions are not renormalizable (or are irrelevant operators) in space- 
time dimensions greater than two. This can be seen by simple power-counting. Con- 
sider the Euclidean action in d dimensions 

S = J tf*z(^-~(^) 2 ) , (1.1) 

where <9 = 7 • <9 and G is the (bare) coupling constant. The coefficient of the four- 
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fermion coupling must have engineering dimension 2 — d: 



G = g I /A d ~ I . (1.2) 

Here, we have defined the dimensionless coupling constant g 2 and have included 
powers of the ultraviolet cutoff A to give the vertex its correct dimension. In d = 2 
the coupling constant G is dimensionless and the model is renormalizable in the 



conventional sense. It was shown by Gross and Neveu [II] that it is asymptotically 
free, has nontrivial scaling and contains fermion-anti-fermion bound states. 

In d > 2, the four-fermion interaction is irrelevant in weak coupling perturbation 
theory. Because of the powers of cutoff 1/A d ~ 2 in the bare coupling, perturbative 
contributions to the full four-fermion vertex vanish as the cutoff is sent to infinity. For 
example, the contribution of each of the bubble diagrams in Fig. 1 (the corresponding 
Feynman integral has ultraviolet divergence ~ A d ~ 2 ) goes to zero like 1/A d ~ 2 if the 
external momentum is kept fixed at a finite value and A — > oo. Thus, the four- 
fermion operator is irrelevant at weak coupling: its contributions to four-fermion 
scattering vanish and its contribution to relevant vertices can be absorbed by local 
renormalizable counterterms. 

On the other hand, it was conjectured long ago |I2| that the four-fermion in- 
teraction could perhaps be relevant if, instead of a simple perturbative expansion, 
one considered some nontrivial resummation of infinite numbers of diagrams. For 
example, if we consider bubble graphs of the kind depicted in Fig. 1, the term with 
n bubbles is 

n-bubble graph ~ g ^+i) A (n+i)(2^d) ^ ^ 

where T(A,p) ~ A d ~ 2 + ... is the contribution from a single bubble. The geometric 
sum over all bubble graphs of the kind depicted in Fig. 1 has the form 

g 2 

sum ~ , , -— - — r + crossed . (1.4) 

A d - 2 - g 2 T(A,p) v ; 

The interaction can be non-zero if 

A d ~ 2 

g 2 = r=— — r + O ((finite mass scales/A) d ~ 2 ) . (1.5) 
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Thus, in order to produce non-zero interactions, the coupling constant must be pre- 
cisely tuned to this strong coupling (ultraviolet stable) fixed point, within errors 
proportional to (finite mass scales/A) d 2 . At the critical point, (|l-4j) behaves like 
sum ~ l/p d ~ 2 so that the effective scale dimension of the composite operator ipip is 1 
(in momentum units) for any d. When d 7^ 2 this differs from its classical dimension, 
(d — 1). Moreover, the dimension of the four-fermion vertex is approximately two, 
rather than its classical dimension 2(d — 1). It is this large anomalous dimension of 
composite operators which makes the four-fermion interaction relevant. 

The latter idea appeared in the work of Wilson |13| who observed that the sum 



of bubble diagrams in (.L4) is the leading order in the large N expansion where N 
is the number of fermion species. He argued that, in dimensions 2 < d < 4, as well 
as the infrared stable fixed point at zero coupling, the four-fermion coupling has an 
ultraviolet stable fixed point at a non-zero value of the coupling constant g 2 > 0. As 
d — > 2 this fixed point moves to the origin, g 2 —* and gives the asymptotic freedom 
of the Gross- Neveu model |TTJ. He also observed that at d — 4 the four-fermion 
theory is trivial and is described by A" non-interacting spinor fields and a free scalar 
field. 

This led to the conjecture that, in 2 < d < 4 the systematic large N expansion 
of four-fermion theory is renormalizable in the conventional sense that, in order to 
remove the cutoff dependence of physical quantities as they are computed order by 
order in 1/N, it is necessary to add local counterterms which have the same form 
as these already occurring in the action [H|||15||[]r5[]. That this is indeed the case in 
d = 3 has only recently been proven using the renormalization group technique of 
constructive field theory [|17] . 



There have also been numerous recent efforts, both analytical |18| and numerical 



191, to understand the critical behavior of the four-fermion model around its fixed 



points. The numerical data in d=3 is in good agreement with analytical calculations 
in the large A" expansion. The ultraviolet fixed point has the interesting feature that 
the critical indices there are much different from those at the Gaussian infrared fixed 
point. 

In order to implement the large N expansion it is convenient to introduce an 
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auxiliary scalar field, so that the action takes the form 

1 , -r , 1 



S = I d a x \ + ~i=UiA + ^0 2 J , (1-6) 



where % = 1, . . . , N. The four-fermion interaction in ( |1.1| ) is recovered by solving the 
equation of motion for or, equivalently, performing the gaussian functional integral 
over (j) in the partition function 

dif>dijjd(j)exp(—S) . (1.7) 

The introduction of results in a Yukawa interaction which, as we shall see, is essential 
to the method we shall introduce in this paper. In the past this trick has been used 
to discuss |20|] the equivalence of the four-Fermi theory and a fermion-scalar field 



theory coupled through Yukawa interactions when the wavefunction renormalization 
constant for the field vanishes. This is known as the compositeness condition 
and, as has been suggested in some recent literature ||10|| , can be used to reduce the 
number of parameters in the fermion-scalar field theory, thereby giving the model 
more predictive power. Indeed, one could consider a generalization of ( |1.6j ) which 
contains a kinetic term for the scalar field \<fi^<fi and possibly higher polynomials - in 
fact all of the relevant operators of the scalar fields theory. Then if one considers the 
renormalized kinetic term, |^0^D0r the compositeness condition can be achieved if 
as a function of g 2 and other parameters, the equation Z<p(g 2 , . . .) = has a solution. 
If, at this solution, correlators of the field theory remain finite, the interacting fermion 
theory which would result from eliminating the scalar fields is also finite. (If the scalar 
field action has terms of higher than quadratic order, the resulting fermion theory is 
rather complicated, possibly non-polynomial.) 

The action ( |1.6| ) is invariant under the discrete chiral transformation 

tpi -> 75lpi , i>i -> -^i75 , -> ~<t> • (1-8) 

where we consider the appropriate definition of 7 5 in less than 4 dimensions. (We 
use spinors with 2^ 2 ' components, where [d/2] is the smallest integer greater than 
or equal to d/2. In any integer dimensions such spinors have an analog of 7 5 which 
is the product of the gamma matrices, n M 7 M -) This chiral symmetry prevents the 



fermion from acquiring a mass. When the four-fermion interaction is attractive and 
strong enough, there is a phase transition and the fermion acquires mass by spon- 
taneous chiral symmetry breaking. The corresponding critical coupling constant is 
just the ultraviolet fixed point. A natural order parameter for the spontaneous chiral 
symmetry breaking is the expectation value of the auxiliary field ((f)) which, in our 
convention, has the dimension of mass. The critical point where the phase transition 
occurs can be approached from either the symmetric or the broken phase. Assuming 
the phase transition here is of second order, as strongly suggested by numerical Monte- 
Carlo simulation [19|, the critical indices calculated in one phase must be equivalent 



to those calculated in the other. 

In the present paper, we shall use the conformal bootstrap to examine four-fermion 
interactions in 2 < d < 4 at the critical point. The conformal bootstrap was invented 
by Polyakov 0, Migdal 0, Parisi || and many others in the early seventies. 

The basic idea is as follows. As the transition point is approached, various quan- 
tities either vanish (such as masses) or diverge and long-range correlations appear. 
Moreover, just at the critical point, the correlation function of fluctuations takes the 
form 



(6s(x)6s(0)) l ^~ ° \x\- {d -' 2 - v) . (1.9) 

This is exactly the form of two-point correlation functions which is implied by scale 
invariance, where the effective scale dimension of the operator 5s(x) is d ~1~ v ■ 



In fact, as was first observed in ||21|| , in addition to scale invariance, at the phase 
transition point the system is invariant under the full group of special conformal 
transformations. To see this, consider the conformal current K®(g 2 ,x) where g 2 is 
the effective coupling constant. Then the divergence of the conformal current is 
related to the beta function of g 2 through the trace of the energy-momentum tensor: 

d tl KZ{g 2 ,x)<xx a m (g 2 ,x) , (1.10) 

which in turn is the divergence of the dilatation current D^(g 2 ,x) and, therefore, 
proportional to the beta function: 

WsV) =d^D„{g\x) oc/% 2 ) . (1.11) 
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As usual, the last equation is understood in the sense of matrix elements. If more 
than one coupling constant is involved, the right hand side of Eq. would be a 

linear combination of all the corresponding beta functions. 

As well as the two-point function, which is actually fixed (up to the critical ex- 
ponent) by scale invariance, conformal symmetry also fixes the form of three-point 
functions up to multiplicative constants. The multiplicative constants are the fixed 
point values of the renormalized coupling constants — by definition — the zero point 
of their beta function. 

To determine the unknown exponents and multiplicative constant (fixed point), 
one invokes the Schwinger-Dyson equations in bootstrap form - the homogeneous 
Schwinger-Dyson equations. The inhomogeneous term arising from bare couplings in 
the Schwinger-Dyson equation for the three-point vertex vanishes since it is multiplied 
by the inverse of an infinite renormalization constant. Since the bare vertex does not 
have the correct form to satisfy the conformal ansatz, this cancellation is in fact 
necessary for consistency of the assumption of conformal symmetry. This is discussed 
in detail in 0. The bootstrap equations for the three- vertex are closed and have a 
conformal solution. They are free of ultraviolet and infrared divergences, as a result 
of renormalization and are sufficient to determine the scale dimension of the fields and 
the value of the coupling constant at the fixed point. In the present model, these are 
the critical coupling constant and the scaling dimensions of the fields ip and <fi. If the 
hyperscaling hypothesis [22] is correct, all other critical exponents can be determined. 

This paper is organized as follows. In Section 2 we analyze the 1/N expansion of 
four-fermion theory at d — 3 to order 1/N, demonstrate renormalizability and review 
the calculation of the anomalous dimensions of ip an d 4> to that order. In Section 3 
we assume conformal invariance of the four-fermion theory at fixed point to derive 
the set of bootstrap equations which determines the critical indices. In Section 4 we 
consider the three-vertex approximation which is legitimate when the index of the 
Yukawa vertex is small, calculate the critical indices and show that results agree with 
the 1/N expansion. In the Appendix we present a derivation of the Parisi procedure 
for the case of the small index of Yukawa vertex. 
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2 1/N expansion 

In this section, we conduct a perturbative expansion of the theory in 1/N with em- 
phasis on the conformal structure. For simplicity, we restrict our computations to 
d = 3, the generalization to any dimensionality in 2 < d < 4 is straightforward. 

2.1 Bubble graphs 

Our convention is that each fermion has two components and the Dirac matrices 
are Hermitean, namely, 7^7^ + 7i/7 M = 25^1 and Tr 1 = 2^ d ^ 2 ^ = 2. In the symmetric 
phase, fermions have no dynamical mass. The Feynman rules in the momentum space 
use the fermion propagator and three-point vertex 

5 (p) = 4 and r = --L, (2.1) 

%V VN 

respectively. 

The scalar, 0, is not a dynamical but an auxiliary field. However, there formally 
exists a bare 'propagator', G, corresponding to the last term in the action ( |1.6| ). G, 
the only free parameter of the theory, has dimension of mass -1 . 

On the other hand, at a fixed point all dimensionful parameters must cancel 
from expressions for correlation functions. Here, the cancellation of the dimensional 
coupling constant G occurs through the the regularization and renormalization pro- 
cedure. If the theory is regularized by introducing a cutoff A which has the dimension 
of mass, one can choose 1/G so that it exactly cancels the terms in the scalar field 
self-energy that are linearly divergent, ~ A. Or, if the theory is regularized by a 
method such as dimensional regularization where linear divergences are absent, 1/G 
set to zero. We can regard both of these as fine tuning of the bare coupling constant 
1/G to the ultraviolet stable fixed point. 

In this Section, we shall introduce a naive cutoff A as a regulator. In this case, 
we can define a dimensionless coupling constant g according to Eq. ( fL.2[) : 

q 2 

G= 9 -. (2.2) 

(In this Section, d=3.) 
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The bare 0-propagator is a momentum independent and is proportional to the 
coupling constant, G. In the 1/N expansion, the fermion bubble chains (Fig. 2) are 
of the same 0-th order as the bare <p propagator. Therefore, the correct propagator 
in the Feynman rules should be an effective one that sums over all one-fermion-loop 
chains. It is referred to as the dressed (ft propagator and its inverse takes the form 

a i/ an 1 f A d 3 k _ k(k + p) A A 1, . . , 

Setting g = g c so that to cancel the linear in A term in Eq. ( |2.3|) , the hi / G—A/ir 2 = 
0, we have the critical coupling constant at this order 



g c = n, (2.4) 



and the dressed propagator 



Ao(p) = ^ ■ (2.5) 

Now we see that, at the order 1/N°, the two- and three- vertices given by Eqs. fl2.lp 
and ( |2.5p have indeed a conformal structure, as if the fields ip and cf) both have 
dimension 1 and the Yukawa interaction is dimensionless. While this value of the 
dimension of the fermion field ip coincides with the canonical one at d = 3, in zeroth 
order in 1/N the scale dimension pg <fi is 1, different from the value 3/2 which is 
implied by the bare propagator .0. This is related to the shift of the scale dimension 
of the composite operator ipip to 1 after summing the bubble graphs and tuning the 
bare coupling to the critical value ( |2.4|) which was discussed in Section 1. Notice that 
the scale dimension 1 coincides with the standard mass dimension of the field 0, so 
that there is no need of a mass parameter in Eq. (|2.5T ). 

2.2 Primary divergences at 0(1/N) 

At the order 0(1/N), we need to consider the diagrams in Figs. 3-6. We shall see 
that the primary divergences are these in the fermion self-energy and in the Yukawa 



1 Since, for large momentum \p\, the momentum-space propagator of a field with the scale dimen- 
sion I is ~ [p 2 ] l ~ d t 2 in a d-dimensional space, the canonical scale dimension of the auxiliary field is 
d/2 while that of a conventional dynamical scalar field is d/2 — 1. 
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vertex. A straightforward calculation of Fig. 3a gives the regulated fermion two-point 
vertex 

-Si f A d 3 k k + p 



! -8i f* d°k k + p 

S (P ' A) = IP -WJ (2nr\k\(k + p) 



= ^ 1 + ^iv ln (f) + finite ]' ( 2 - 6 ) 

and of Fig. 4a the regulated three-point vertex 

1 8 f A d 3 k {k + pi){k+p 2 ) 



N VIPJ (2ir) 3 \k\(k + Pl ) 2 (k + p 2 ) 2 
' t 1 - Jii? finite], (2.7) 



v^V L tt^N K pl 

where p^ ax is the largest of p\ and p|, with and p 2 the incoming and outgoing 
fermion momenta, respectively. 

The diagrams in Fig. 5, which have a subdiagram with a fermion loop and three 
boson legs. The subdiagravanishes because of parity symmetry. We shall see in the 
next section that this form of the vertex function is in fact prescribed by conformal 
invariance. 

To render the above two- and three-point vertices finite, we need to introduce 
local counterterms, Figs. 3a' and 4a', 

Fig.Sa' = -^ln(^)^, (2.8) 

Fig.4a' = -— l=ln(^)#V, (2.9) 

respectively, where a reference mass scale \x has been introduced. The corresponding 
renormalization constants, using minimal subtraction, are 

Z *< A > = l ~^N^f' < 2 - 10 > 
The anomalous dimension of the fermion field is 
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From (|2.7|) , we see that the three-point vertex, while dimensionless at the tree level, 
develops an anomalous dimension, —27, beyond the tree level. We shall call 7 the 
index of the Yukawa vertex. At the order 1/N, 

7=-^- ("3) 
2.3 A demonstration of renormalizability 

The renormalization of the operator 2 is crucial. By power counting, the inverse 
of the two-point function (00) has engineering dimension one, instead of two, which 
would be the case for a dynamical scalar field. The self energy for the scalar is linearly 
divergent, with the divergent terms ocurring in Feynman diagrams cancelled order by 
order in the 1/N expansion by adjusting the value of the critical coupling constant. 
The remainder of the scalar self-energy can (and generally does) contain individ- 
ual integrals with non-local divergent terms like \p\ ln(A 2 /p 2 ). It should be rendered 
finite by the insertion of the fermion wave-function renormalization and vertex renor- 
malization counterterms into the internal lines of lower order diagrams. Otherwise, 
non-local counterterms would be necessary to cancel ultraviolet divergences and the 
1/N expansion of the theory would be non-renormalizable in the conventional sense. 
In this Section, we shall demonstrate that, to order 1/N, these nonlocal divergences 
are indeed cancelled by local counterterms. 

To see this, we consider the diagrams at the order 0{1/N), Fig. 6. Figs. 6a and 
6b each has a subdiagram, the one-loop fermion self-energy, that is logarithmically 
divergent, as calculated above. Therefore accompanying Figs. 6a and 6b are Fig. 6a' 
and 6b' containing the counterterm fl2.8|) . Similarly, Fig. 6c has a subdiagram, the 
one- loop vertex correction, that is also logarithmically divergent, and therefore we 
need to calculate Figs. 6c' and 6c", together with Fig. 6c. Once these counterterm 
diagrams have been taken into account all logarithmically divergent terms cancel 
out between Figs. 6a and 6a'; Figs. 6b and 66'; and among Figs. 6c, 6c' and 6c". 
Consequently, no momentum dependent counterterms are necessary. 

The linearly divergent terms in the <fi self-energy can be cancelled by adjusting 
the critical coupling constant g 2 , which has been determined to leading order in 
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the previous subsection. By this analysis, we have confirmed the renormalizability of 
( |1.6|) at the order 1/N. In the language of multiplicative renormalization, introducing 
counterterms implies introducing renormalization constants for the vertices appearing 
in the action. In the the present case, for is not independent but takes the form 
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zt = ! 1 + 3^v"V]- (2 - 14 > 



The following explicit calculation confirms ( |2.14| ): 



Fig. 5a + Fig. 56 + Fig.5c 



8 f d 3 k d 3 q k(k + p)k(k + q) k(k + p)(k + p + q)(k + q) 



NJ (27r) 3 (27r) 3 l k\k + p) 2 (k + q) 2 q k 2 (k + p) 2 (k + q) 2 (k + p + q) 2 q> 
2 f A d 3 q r 2 q ■ p p , 



NJ (27r) 3 \q + p\q q 2 \q + p\p Q ■ p\Q + p\ 



2 



1 r . 2 , A , 

[- 2A + -pln(— ) + const x p] . (2.15) 



tx 2 N v r V 

By using (|2.3|), we have 



A- 1 ^, A) = A - [A _ h _ J_ [ _ 2A + |pln(^) + const x p] . (2.16) 

Choosing g = g c so that A/g 2 — (1 — 2/N)A/it 2 = 0, we have the critical coupling 
constant at the order 1/iV 

^ = vr(l + A).^ = 1(1- A). (2.17) 

A* 7T A 

The renormalized scalar two-point function is 

A- 1 ^) = Z,(A)A-\p, A) = |[1 - AA-ln(^)] , (2.18) 



where we have used the wavefunction renormalization constant Q2.14|) . The anoma- 
lous dimension of is 



From (|2.12 ), ( |2 . 1 3|) and (|2.19|) , it can be checked that the index of Yukawa in- 



teraction and the anomalous dimensions of if) and at the order 0(1/N) satisfy a 
relation 
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7 = 7^ + \l4> ■ (2-20) 

As we shall show in the next section, this is not an accident but a result of conformal 
symmetry. The point is that the renormalization constants in the conformal theory 
are equal to their fixed point values: 

zm = {^y\ (2.2i) 

*i(A) = (^)\ (2.22) 

Z,(A) = (^) 70 - (2.23) 

The expressions (|2 . 10|) , (|2.11|) and Q2.14| ) which are explicitly calculated at the order 
0(1/N) can be viewed therefore as the 0(1/N) terms of the expansion of (|2.21|) , 
(|2T22|) and ( ^23|) in 1/N with ^, 7 and 7^ given by Eqs. ( gl2|) , (|2T3|) and (|2~T9l) , 



respectively. 

3 Conformal bootstrap 

In this section, we derive the bootstrap equations for the system with four-fermion 
interactions in d dimensions. The result is three equations which determine two 
anomalous dimensions and the value of the critical coupling constant. We start with 
a systematic analysis of the conformal structure of a general theory involving fermions 
and bosons with Yukawa interaction (s). 

3.1 Conformal structure 

Let I and b be the scaling dimensions of the primary spinor and scalar fields, respec- 
tively. Under a scaling transformation of the coordinates, x — > px, the transformation 
of the fields are defined as 

ijj'(px) = p~ l ijj(x) , ip'(px) = p~ l ^j{x) , 4>'{px) = p~ b (j){x) . (3.1) 
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Moreover, the special conformal transformation of the fields are defined as: 

ijj{x) -> ijj\x') =a l - 1/2 {l + tx)tfj{x) , (3.2) 
${x) -> = i - 1 / 2 $(x)(l + xf) , (3.3) 

_> <p\x>) = a b x <P(x), (3.4) 

where <r x = cr t (x) = Ij^l 1 ^ = l+2t-x+t 2 x 2 , under a special conformal transformation 
of coordinates, parameterized with a constant vector t M , 

X " X " _ i + 2t -x + tW ■ ^ 

Scaling and special conformal transformations together with translations and Lorentz 
rotations comprise the conformal group which has (d + l)(d + 2)/2 generators in d 
dimensions. 

The transformations of the primary fields determine the transformations of the 
Green functions. For instance, under the conformal transformation Q3.5|) , the two- 
and three-point functions transform as 

G'(x' - y') = W)> = (^^)^ 1/2 (1 + tx)G(x -y)(l + yt) , (3.6) 

D\x' - y') = (7>V)W)> = (cr x a y ) b D(x - y) , (3.7) 
G' 3 (x',y';z') = {T<t>>{z'W{x'W{y')) = a b z (a x a y ) l -^ 2 (l + tx)G 3 (x, y;z)(l + yt) .(3.8) 

The unique (modulo normalization) solution to the special conformal transformations 

G{x ~ y) = a^-WN(1) [(x - y) 2 } l+1 l 2 ' (3 ' 9) 
D (?~v) = 4 »-<WV(6) \(x-yf] b ' (3 ' 10) 
G 3 (x,y;z) = C ^ _ z ^ b/2+1/2 ^ _ y)2 ]f/2+i/2 ^ _ y )2]j-6/2 ' ( 3 ' U ) 



where /z = d/2, and 



v ; r(r + l/2) ' v ; T(r) v 7 
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Note that, unlike the case of d=2 where the conformal group has an infinite number 
of generators, the finite dimensional conformal symmetry in d > 2 restricts but does 
not fix the functional forms of higher point correlation functions. 

( |3.9| ) and ( |3.10| ) are normalized so that their Fourier transformations take a simple 
form 

G(p) = -Ap 2 ) l ~ h+1/ \ (3.13) 
ip 

D{p) = ^(p 2 ) b - h+1 . (3.14) 



Amputating the external legs i from the three-point function (|3.11 ), we obtain the 
conformal three-point vertex (depicted in Fig. 7) 



T(x,y;z) = g^-N {l )N\ h / 2 )N{l - b/2) 

x — z z — y 1 

' {(X - Z )2]h-b/2+l/2 [7^ _ y )2]h-6/2+l/2 [7^ _ y^h-l+b/2 

which is normalized so that the Fourier transformation of ( |3.15| ) is 



(3.15) 



f d d k k + pi k + p2 1 

T(p 1 ,p 2 )= g ,T(h)N( 1 )j ^ [(fc+pi)2f/2+1/2 [(Jfe + P2 )2]V2 + i/2 rrbt 2 . (3- 16 ) 

where g* is the unknown dimensionless fixed-point coupling constant and 

j = l + b/2-h , (3.17) 

is defined as the index of the vertex. By dimensional analysis, it is easy to check that 

the dimension of the vertex is —2 / ~f. 

In the previous section, we have introduced in d = 3 the anomalous dimensions 
and 7^ by Eqs. ( |2.21| ) and ( |2.23| ). They are related to the scale dimensions I and 

b as follows: 

/ = h-l + j^, (3.18) 
b = 1 + 70- (3.19) 



2 This amputation in conformal theories replaces the dimensions by 'shadow' dimensions I* = d—l, 
b* =d-b. 
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When these formulas are substituted into Eq. (|3.17 ), it recovers Eq. Q2.20Q which is 



obtained in the previous section by explicit calculations at the order 0(1/ N). It is 
worth also noting that when ( 3.18|) and ( |3.19| ) with 7^, = 7^ = are substituted into 



( |3.16| ) the result coincides in d = 3 with the vertex function ( |2.7j ) which explains its 
conformal origin. 

3.2 Bootstrap equations 

The critical coupling and the indices b and / are determined by the set of homoge- 
neous Schwinger-Dyson equations, which we shall now derive. To simplify notations, 
we introduce operators ST and ST as 

( Pl \ST\p 2 ) = G( Pl )T( Pl , P2 )G(p 2 ) , ( Pl \ST\p 2 ) = G( Pl )T( Pl , P2 )D( P2 - Pl ) , (3.20) 

where the matrix multiplication over the spinor indices is implied. Here, S and 5" 
are the operation of attatching two fermion or a fermion and boson leg to the vertex 
function, respectively. Then the set of Schwinger-Dyson equations reads 

T = go + TSK, (3.21) 
E^G^-G- 1 = g S'T, (3.22) 
U^D^-D- 1 = -g NTr(ST), (3.23) 

where Tr is the trace over the spinor indices and K stands for the standard Bethe- 
Salpeter kernel. 

The bare coupling g in Eq. ( |3.21| ) must vanish so that this equation is consistent 
with the conformal ansatz. This is realized by renormalization 0. Moreover, since 
the three vertex is fixed by conformal invariance, both sides of ( |3.21| ) must have 
the same dependence on the space-time coordinates. Therefore the vertex bootstrap 
equation takes a simple form 

l=g 2 J(l,l,b;g«), (3.24) 

where /, named the vertex function, is a function solely of the scaling dimensions 
carried by the legs and of the critical coupling. The vertex function, to be explored to 
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some extent in the next section, plays a key role in the conformal bootstrap approach, 
calculated below, lines is shown I or b). 

However, one cannot immediately set the bare coupling go to zero in Eqs. ( |3.22|) 
and ( |3.23| ). The reason is that the quantities that go multiplies in ( |3.22| ) and ( |3.23| ) 
are divergent. Therefore the right hand side of ( |3.22| ) and (|3.23| ) has the form ■ oo. 

To deal with this uncertainty, we differentiate both sides of Eqs. ( |3.22j ) and ( |3.23| ) 
with respect to the external momentum k. Denoting Jp- = S Q and Jp- = Tl a and 
replacing go with V — TSK and V — TS'K', we have 

s Q = (r - rs'K')s' a r + (r - rs'K')s'r a , (3.25) 

U a = -NTr[(T-TSK)S a T + (T-TSK)ST a ] , (3.26) 

where K' denotes the boson-fermion Bethe-Salpeter kernel. The skeleton expansions 
of K and K' are depicted in Fig. 8. The vertex bootstrap equation now is 

T = TSK = TS'K' . (3.27) 



Now we see that the first term on the right hand side of Eqs. (|3.25| ) and (|3.26|) 



vanishes as the bracket is zero, according to fl3.27|) , and S' a T and S a T are finite. 
However, the second term of Eqs. ( |3.25|) and (|3.26|) is still the form ■ oo. 



3.3 Parisi approach 

To remove the remaining uncertainty in ( |3.25| ) and ( |3.26| ), we use a regularization 
proposed by Parisi 0. The idea is to shift slightly the dimension of one of the 
operators in the Green functions: 

l^l' = l + e , or b -> b' = b + e . (3.28) 

The corresponding change of the expression on the right hand side of the vertex 
bootstrap equation is 

rS'K' = g 2 J(l + 6,l,b ig *)r , (3.29) 

rSK = g 2 J(l,l,b + e-g*)r , (3.30) 
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where / is defined similarly to the right hand side of Eq. (|3.24 ). Using Eq. (|3.24 ) we 
obtain 

2 df(l',l,b;g*) 



- 9* 



dl' 



ers' r r 



IE 



-Ng 



2 df(l,l,b';g* 



db' 



eTr [rST a } . 



(3.31) 
(3.32) 



b'=b 



Notice that, in this procedure, the dimension of only one external line is analytically 
continued in the equation. 

The calculations of r e S"r o and r e ST a are performed similarly to those of the 
scalar $ 3 theory || and of $ 4 theory [Q. In the fermion model, a useful formula is 
Mi 

dz (x\ — z)(z — x 2 ) 



71 [( Xl - z) 2 ] Si+ 2 \(x 2 - z) 2 } S2+ 2 [(x 3 - zyy* 

x 13 x 23 N(8 l )N(8 2 )N(5 3 ) 



\x 2 13 ] h - 52+1 2 \x 2 ^ h - 5l + ^ 



-23 J 



[x 2 l2 ] h ' S3 



(3.33) 



It is valid providing 

Si + 8 2 + 5 3 = d . (3.34) 

The meaning of these formulas is very simple: if one performs the conformal transfor- 
mation ( p.5|) of x's, the left hand side of Eq. ( |3.33|) obeys Eq. (|3.8|) providing Eq. (|3.34j) 
for 5's is fulfilled. Therefore, the integral is determined (modulo an ^-independent 
factor) by the conformal invariance to be the right hand side of Eq. (|3.33| ). 

Using Eq. (|3.33| ) in calculating r e 5"r Q and r e 5T a for each of the integrals (the all 
are of the type ( |3.33|) ) and picking up the terms ~ e -1 , we rewrite Eqs. ( p.31|) and 
(p|) as § 



A^Trl 



(4tt) 

at 

(47T)' 



-Y{h)N{ 1 )N{l)N{d - l)N 2 {b/2)N{l - 6/2) 



df(l',l,b;g* 



■T(h)N(-f)N(b)N(d - b)N(l - b/2)N 2 {b/2) 



dl>/2 
df(l,l,V;g t 



,(3.35) 



i'=l 



dV/2 



b'=b 



(3.36) 



Eqs. ( p.35|) and ( p.36| ) together with Eq. (|3.24|) determine the critical indices. All the 
three equations depend on the vertex function f(l, I, b), which we shall discuss in the 
next section. 
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4 Three-point vertex function 



The skeleton expansion of the Bethe-Salpeter kernel K (or K') in the vertex boot- 
strap equation ( |3.27| ) contains an infinite number of terms. If the critical coupling 
is reasonably small, the three vertex approximation makes sense. Then the vertex 
equation has a solution with small g* for 7<1. The point is that the integral in the 
vertex equation is slowly convergent if 7 <C 1 and is therefore large (~ 1/7). This 
enables it to cancel the small factor of g 2 in the bootstrap equation. 

4.1 Momentum-space analysis 

It is convenient to perform calculation in the momentum space where the term ~ I/7 
comes from the region with large integration momentum. The Fourier image of the 
conformal vertex ( |3.15| ) has been given in ( |3.16| ). The integration on the right hand 



side of ( |3.16| ) has an expression [2~5] of the Appel functions F4 in general, and takes a 



simple form for 7 <C 1. We shall derive the corresponding formula directly analyzing 



the integral on the right hand side of Eq. (|3.16|) 



Let us first note that the coefficient in front of the integral is ^(7) « 7 for small 
7, so that one has to peak up the term ~ I/7 for the vertex to be of order 1. It is easy 
to see that this term comes from the region of integration with \k\~> max{|pi|, \p2\\- 
(Recall that \p\ — p 2 \ < max{|pi|, |p 2 |} in an Euclidean domain.) One gets 



d d k k + pi k + p 2 1 

\{k + Pl ) 2 ] 6 /2+l/2 p + p2 )2]6/2+l/2 [ fc 2]J-6/2 

dk 2 1 



c{ P \,pl} [k 2 ] 1+ i 7 (max{plp 2 2 }y 



(4.1) 



and 

T{p ^ ) = 9 % m a X {pip 2 }y - (42) 

This dependence of the three-point vertex solely on the largest momentum is 
typical for logarithmic theories in the ultraviolet region where one can put, say, p\ = 
without changing the integral with logarithmic accuracy. This is valid, however, if 
the integral is fast convergent in infrared regions. For our case this means that the 
dimensions / and b should be far away from the values at which the integral on the 
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left hand side of Eq. ( |4.1|) is infrared divergent. For instance if I — 6/2 m h, the term 
of the form 

infrared term = [ _ b/2 _ h [ max {p2 )p 2}p/2+i/2[ mm { p 2 ;p 2} ]i+ i/2-fc > ( 4 ' 3 ) 

were appear in the integral which would depend therefore both on the maximal and 
on the minimal momenta. 

In particular, for 7 = -2/(Ntt 2 ) Eq. flO|) reproduces to order <D(1/N) Eq.(2.??) 
which was obtained by direct calculations. For this reason our analysis of the case of 
small 7 is very similar to that of Section 2. 

Let us explicitly write down the dependence of the vertex function / on g*: 

f(l, I, b; g.) = g 2 J 3 (l, I, b) + g 2 J 5 (l, /, b) + . . . , (4.4) 
where f%(l,l,b) is associated with the three-point vertex contribution on the right 



hand side of the vertex bootstrap equation (|3.27|) , f 5 (l,l,b) with the five- vertex terms 
and so on. If one restricts oneself to the three-point vertex term, the resulting equation 
is depicted in Fig. 9. This truncation is referred as the three- vertex approximation. 

It is not difficult to calculate / 3 if 7 <C 1 . In order to pick up the term ~ 1 /j one 
performs quite similar calculations to the logarithmic case. Let p\ > p\. Then the 
three vertex equation reads 



1 2 f d d k 1 k+pi 

- -9* 



[p|]7 y *7 (2n) d [max{pl k 2 }p [(k + Pl )^h-i+i/2 

1 k + p 2 1 1 



(4.5) 



[max{pl k 2 }p [(k + p 2 y]h-i+i/2 [ max {pl k 2 }p [{k 2 )] h - b 

where the momenta are depicted in Fig. 9. 

The term ~ I/7 on the right hand side of Eq. ( |4.5D comes from the domain where 
the integration momentum k 2 > p\. Then we have finally 

I ^00 dk 2 I 
f 3 (l, l,b)=- ^ hT{K) j pl p|TT7 = " (A<K) h T{h) 7 ' (4 ' 6) 



It has been verified simultaneously that the momentum dependence (O) is repro- 



duced by the three vertex equation at small 7. This is a consequence of conformal 
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invariance which fixes as is mentioned above the coordinate (or momentum) depen- 
dence of the three-point vertex. 

To apply the procedure described in the previous section, we repeat the calculation 
above with the dimension of one of the external lines shifted and have 
x m i n 1 , I' + l + b-d 

A(i ''- 6) = ->= 2 ; (47) 

„ /, , ,/n 1 ,, I + I + b' — d , 
MM = 7= 2 ' (4 ' 8 > 

Such a modification is prescribed by the dimensional analysis. 

Substituting Eqs. (|4.7|) and (|4.8|) into Eqs. ( |3.35|) and ( |3.36|) for propagator, we 



arrive at the algebraic equation set that determines the critical indices in the three 
vertex approximation: 

1 = ^* (49) 

(47r)T(/i) 7 ' 1 ; 

1 = ^K^(b/2)N(l-b/2)N(d-l), (4.10) 
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2 



-iVTrl-^— N 2 (b/2)N(d-b) . (4.11) 

(471")" 



4.2 The solution at large N 



We look for a solution to the algebraic equations (|4.9|) - ( |4.11|) for small 7 = l+b/2 — h. 



Let us assume the scaling dimension of ip is close to the canonical value, I ~ h— 1/2, 
i.e. 7^, which is defined by Eq. Q3.18 ) is small. In order that 7 = 7^ + (6 — l)/2 <C 1, 



the scaling dimension of </> should be b ~ 1, i.e. 7^ which is defined by Eq. (|3.19|) 
should also be small. 

From Eqs. (O) and (|4.11|) we have Let us first equal the r.h.s.'s of Eqs. (|4.11|) 



and (|4.10|) for b ~ 1 and l-/i + l/2<l. Solving the equation, we get 
, , . 1 T(2h - 1) sin \ir(h - 1)1 

'^-'^mVwTIF' (4 ' 12) 



and from Eqs. (Q and ( ^T0|) 



1 r(2fc-l)sin[7r(fe-l)] 

7 " iVTrl r 2 (/i)7r ' { ' 
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which gives 

6-1 2 r(2ft)sin[,r(fe-l)] 

iVTrl r(/i)r(/i + 1)tt ' 1 ' ' 

Now we see that the scaling dimensions of ip and are such that 7^ and 7^ are 
~ 1/iV and therefore are indeed small in the large iV limit. Moreover, since g\ ~ 1/iV, 
it is verified that for large iV the three-vertex approximation is a good one. 
Substituting d = 3 and Tr 1 = 2 into ([OP , (gig) and ( PUD , we obtain 

' = i + jp*. < 4 ' 15 > 

These values of scaling dimensions coincide with those, given by ( lip and (gig) , 
which are calculated in Section 2 by the 1/N expansion while the (nonuniversal) value 
of g*, which is determined by Eq. ( f4.9|) , is to be compared to fl2.17|) . 

Notice that the value of 7^, which is given by the last term on the right hand side 
of Eq. Q4.12 ) is positive for 2 < d < 4. In other words, the dimension I is bigger 
than the canonical value. This is in agreement with the Lehman theorem which says 
that the momentum space propagator in an interacting theory without ghosts should 
grow with momentum faster than the free one. For the auxiliary field the Lehman 
theorem is not applicable since there is no pole term in the spectral representation 
of the propagator, so that the dimension b can be arbitrary in four-fermion theory. 
However, if we were consider instead a field theory with dynamical scalars coupled 
to fermions through Yukawa interaction, it would be described by exactly the same 
bootstrap equations with the solution given by (|4.12|) - ( [4.14 ). In this case the Lehman 
theorem for the scalar field would say that b should be bigger than h — 1, i.e. than 
the canonical value for a dynamical scalar field. We see that Eq. (|4.14j ) obeys this 
restriction providing d < 4. Therefore, our solution is applicable to Yukawa theory 
as well as to four-fermion theory. This confirms that the two theories are equivalent, 
as was discussed in Section 1, for 2 < d < 4. 
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4.3 Estimating high-vertex contributions 



Let us stress once more that the system of equations ( [4.9| ) - ( Oil ) is valid provided 



no infrared divergences emerge in the vertex bootstrap equation. This is indeed the 
case for the values of /, 7 and b given by Eqs. ( |4.12j ) - ( [L14j ) (or ( [4.1 5| ) and ( (4.16p for 
d = 3). 

The three vertex approximation works if five etc. vertex contributions are small. 
It is easy to estimate the order in 1/N of these higher- vertex terms which have been 
disregarded. This can be done as follows. 

The index of integrals in the Bethe-Salpeter kernel (which coincides with that for 
a four-fermion vertex) is 

t = 2l-h. (4.17) 

If t 3> 7, five etc. vertex contributions are suppressed as 

A 2 7 

— ~ 1 . (4.18) 

This is true in particular for d = 4 - e (e < 1) when 7 ~ e but t ~ 1. However, 
if t ~ 7, all parquet graphs are of the same order and should be summed up (like 
four-boson interaction in 4 — e (e < 1). Another example is Gross- Neveu model in 
2 + e (e<C 1) when the s-channel graph is compensated by the w-channel one so that 
the three vertex approximation works while t ~ 7. 

5 Conclusions 

The method of conformal bootstrap turned out to be very useful for calculating the 
critical indices for the four- Fermi interaction at the fixed point. From the viewpoint 
of the 1/N expansion this method corresponds to a resumming of the perturbative 
series, so that at each order in 1/N one deals with manifestly conformal invariant 
two and three point functions whose dependence on momenta is fixed by conformal 
invariance. The only unknown quantities: the scale dimensions of ip and ipip field 
as well as the critical value of the coupling constant are determined order by order 
in 1/N. In d = 3 we have shown by explicit calculations how the logarithms of 
perturbative theory are nicely collected into conformally invariant structures. As is 
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discussed above, the reason behind this is the exponentiation of the log's into scale- 
invariant expressions which is a nonperturbative phenomenon from the viewpoint on 
the 1/iV-expansion. 

We have demonstrated that the calculation of critical indices by the method of 
conformal bootstrap is an economical one in the sense that the number of skeleton 
graphs in the vertex bootstrap equation is much less than the total number of graphs 
of the perturbative expansion to the given order in 1/N. Moreover, one only needs 
to calculate the vertex function while the two propagator equations are expressed via 
it. 

It would be very interesting to extend the conformal bootstrap approach to the 
current-current four-Fermi interaction and to quantum electrodynamics. In these 
cases one might expect further simplifications due to the conservation of the vector 
current. 
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Appendix A Parisi procedure at small 7 



Here we derive the conformal boostrap equations from Schwinger-Dyson equations 
for the self-enegies, for small 7. According to ( |3.31| ) and ( |3.32| ), we should extract 
the e _1 term in r e ST a and r e S"r Q . To the leading order in 7 we have 

dq 



Tr 

dk 



dq ghT{h) 



Tr[P5rj =- Tr j ^r(q,p + q)G(p + q)-^-T(p + q,q)G(q) 



q +p 



7v h max{q 2 jP 2 }^/ 2 [qtf-i+W [( g + p )2]fc-i+i/a 
1 k d k — p 



[{q - k) 2 ] l - b / 2 [fc 2 ] fe /2+i/2 g pa _ p )2]6/a+i/2 



Since the integral over dq equals 2ii h /eT(h), it turns out to be 



(A.l) 



Tr [T e ST c 



2 2 d r dk k k — p 

I 9 * 1 T dp~J ^ [A; 2 ] fc /2+l/2 [( fc _ p )2]b/2+l/2 

\gh^^y] h ~ b N{d-b)N\\). 



(A.2) 



Similarly, we obtain 



rs'v a = \gh^Mp 2 ) h - l - ll2 ]N{\)N{d - i)N{i - b -). 



(A.3) 



Substituting Eqs. ( |A.2| ), ( |A.3| ) into ( |3.31| ) and Q3.32| ), we arrive at, in the small 
7-limit, Q, piq ) and ( pip . 
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Figures 




Fig. 1 The four-Fermion scattering in the four-Fermion theories. The solid line 
is the fermion propagator. 




Fig. 2 The summation of infinite series of one-loop Fermion bubble chains gives 
the dressed scalar propagator, which is of order O(N ). The cross line is the 
dressed scalar propagator, the dashed line is the 'bare' scalar propagator. 
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Fig. 3 The fermion self-energy at 0(1/ N). a' is the proper counterterm. 




Fig. 4 The vertex correction at order 0(1/ N). a' is the proper counterterm. 
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Fig. 5 The null diagrams in the symmetric phase at order 0(1/N). 




Fig. 6 The scalar self energy at order 0(1/N). a', b', c' and c" are due to the 
proper counterterms. 
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Fig. 7 The conformal Yukawa vertex constructed by truncating the external legs 
of the conformal three-point function ((pipijj). 
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Fig. 8 The skeleton expansions of a) Fermion-Fermion and b) boson-Fermion 
Bethe-Salpeter kernel. The double solid line is the full Fermion two-point func- 
tion, the double dashed line is the full scalar two-point function. 
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Fig. 9 The three-vertex approximation of the vertex bootstrap equation. 
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